Topological invariants for the fractional quantum Hall states 
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We calculate a topological invariant, whose value would coincide with the Chern number in case 
of integer quantum Hall effect, for fractional quantum Hall states. In case of Abelian fractional 
quantum Hall states, this invariant is shown to be equal to the trace of the Ji'-matrix. In case of 
non- Abelian fractional quantum Hall states, this invariant can be calculated on a case by case basis 
from the conformal field theory describing these states. This invariant can be used, for example, to 
distinguish between different fractional Hall states numerically even though, as a single number, it 
cannot uniquely label distinct states. 

PACS numbers: 67.85.Lm, 03.75.Ss, 67.85. Hj 



It is well known, following the pioneering work of 
Thouless and collaborators [1], that the Hall conductance 
of a free fermion system is an integer valued topological 
invariant, in units of /h. In recent years the concept 
of topological invariants was generalized to a far larger 
variety of free electron systems, termed topological insu- 
lators [2] , and a plethora of topological invariants for free 
fermions in a variety of spatial dimensions was proposed 
and classified [3] . These invariants are usually written in 
a form specifically adapted to noninteracting fermionic 
systems. For example, they are typically expressed in 
terms of single particle Bloch waves of the underlying 
noninteracting Hamiltonians. The topological invariant 
for the integer quantum Hall effect is the Chern number 
characterizing the bands of free fermions moving in a two 
dimensional space. 

At the same time, it is now well understood that these 
topological invariants can be reexpressed in terms of sin- 
gle particle Green's functions [4]. In this form, they are 
defined even if interactions are switched on. Their exis- 
tence reflects the topology of the Green's functions. How- 
ever, a topological invariant written in terms of Green's 
functions no longer corresponds to a response to an ex- 
ternal perturbation. For example, the Chern number 
reexpressed in terms of Green's functions, once the in- 
teractions are turned on, is no longer necessarily equal 
to the Hall conductance (although they remains equal if 
the interactions are weak [5] , within the integer quantum 
Hall state). Nevertheless, they retain certain physical 
meaning thanks to a relationship between the topologi- 
cal invariants in the bulk and at the edge, first derived by 
G. Volovik in case of the Chern number type invariant [6] , 
and subsequently generalized by to a larger class of in- 
variants in a variety of spatial dimensions in Ref. [7] . The 
edge of a topological insulator is gapless and therefore is 
not an insulator; the edge invariant is of the type used 
to characterize topological (semi)-metals such as the one 
studied in Ref. [8]. The precise correspondence between 
the edge and the bulk invariants is described in Ref. [7]. 

Here we will use this bulk-boundary to calculate the 
topological invariant, the Chern number reexpressed in 



terms of single particle Green's functions, in a variety of 
fractional Hall states, where it is by no means equal to 
the Hall conductance. We will see that, quite generally, 
for the states described by a if-matrix [9] , this invariant 
is equal to the trace of that matrix. In a Read-Rezayi 
non- Abelian state the invariant is equal to A/ -I- 2, where 
M is defined as in Ref. [10]. In particular, in the Moore- 
Read (Pfaffian) state [11] at the filling fraction v = 5/2, 
M = 1 and the invariant is equal to 3. Finally, for the 
anti-Pfaffian state at z/ = 5/2 the invariant is simply 1. 

The method we use for calculating the invariant re- 
lies on the detailed knowledge one typically has about 
the low energy theory, and thus the Green's functions, 
of the edge of fractional Hall states. Indeed, the Green's 
functions in the bulk are generally not known, thus eval- 
uating the invariant directly in the bulk does not seem 
to be possible. However, the edge topological invariant 
relies on the knowledge of the edge Green's functions and 
can be evaluated directly, which is what we do here. The 
bulk-boundary correspondence states that the edge and 
bulk invariant are equal, and thus the edge calculation 
directly produces the value of the bulk invariant. 

The utility of this observation, in our opinion, lies pri- 
marily in the possibility of evaluating the invariant nu- 
merically. Indeed, if a quantum Hall state candidate is 
found numerically by exact diagonalization, for example, 
in order to facilitate its identification one may be able 
to calculate the Green's functions, and the invariant, in 
this state. While the invariant does not uniquely label 
different quantum Hall states as follows from its values 
listed above, it provides an additional test which helps 
with the identification of the state found approximately 
numerically. Note that here we evaluate the invariant 
by looking at the edge theory of the corresponding frac- 
tional Hall state. Numerically it would not be convenient 
to look at the edge theory directly, rather it is more con- 
venient to construct the theory on a sphere, for example. 
Then a direct evaluation of the invariant in the bulk be- 
comes possible, which could be compared with the theo- 
retical evaluation via the edge theory. 

This invariant may also prove useful in the search for 
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fractional quantum Hall analogs in interacting topologi- 
cal insulators in three spatial dimensions. 

To proceed, let us state the relationship between the 
bulk invariant and the boundary as it stands in two di- 
mensions. Given the Green's function for an infinite two 
dimensional system Gab{i^,kx,ky), where oj is the Mat- 
subara frequency, kx and ky arc two dimensional mo- 
menta, and indices a, b refer to the bands and/or spin 
and flavor of the fermions, the topological invariant is 
written as 



-^2 = ^ ea/37 
a/37 



dtocPk 
247r2 



tvG-^dk^GG-^dk,GG-^dk,G, 

(1) 

where a, and 7 are summed over 0, 1 and 2, ko = oj, 
ki = kx and k2 = ky. The subscript 2 in A'2 indicates 
that this invariant is for two dimensional systems (suit- 
able generalizations exist in other dimensions). 

This number is an integer regardless of the origin 
of Gab{'^^kx,ky). In the absence of interactions, G = 
[iw — 'H]~, where 'H is the Hamiltonian of noninteract- 
ing fermions. Then N2 is equal to the combined Chern 
number of the negative energy, single particle bands, and 
thus coincides with the Hall conductance. In the pres- 
ence of interactions, it is no longer generally equal to Hall 
conductance, but retains its topological nature (remains 
strictly an integer). 

The bulk-boundary correspondence can be stated as 



N2 = No{A) - No{-A), 



(2) 



where A^o is the edge topological invariant calculated for 
the system with a single edge as 

No{k) = j jdx'dx"tvG-^{w,x',x",k)x 

duiG{uj,x" ,x' ,k), (3) 

where G{lo,x' ,x" ,k) is the Green's function of a system 
with a single edge, which depends on the momentum k 
along the edge and the two coordinates x' and x" perpen- 
dicular to the edge (since the system with an edge is not 
translationally invariant in the perpendicular direction), 
and where A is a suitably chosen (sufficiently large) pa- 
rameter. For the purpose of Eq. (3) the inverse Green's 
function G~^ is defined by 

/ dx'G~^{uj,x,x',k)Gbc{to,x',x",k) = 5acS{x-x"). 
b •' 

(4) 

The edge invariant Nq encodes the topological infor- 
mation about the edge (in fact, signifies that the edge is 
a topological metal [6]). For example, in the absence of 
interactions iVo(A) — No{—A) is just the number of chiral 
edge modes, so that Eq. (2) simply reflects the fact that 
the number of chiral edge modes is equal to the bulk con- 
ductance N2, something which is well known and can be 



established in other ways. In the presence of interactions, 
the meaning of Eq. (2) is a little less transparent and was 
established in Ref. [12] to be the difference between the 
number of chiral edge modes and chiral edge zeros. 

To further clarify this point, observe that if one intro- 
duces the eigenvalues A„ of G, 



^ / dx'Gab{uJ,X,x',k)4"\x')=Xnilj^:'Hx), 



then 



(5) 



(6) 



which is the sum of windings of the phases of A„ as w 
goes from —00 to 00 divided by 27r. This last statement 
can be used to calculate iVo(fc) easily. In particular, in 
the absence of interactions we have A„ = l/{iui — e„(/c)), 
where e„(fc) are the energy levels of the system with an 
edge and with momentum k along the edge, and No{k) 
can be calculated as 



No{k) = ^arg A„|^^ 



c<;=oo 
00 



signcnik). (7) 



Note that A^o does not evaluate to an integer, because of 
the slow decay of A„ with ui. However, diff'erences will 
still take integer values. Thus No{A) — No{—A) counts 
the number of energy levels whose energy changes sign as 
k is varied from —A to -|-A. This is simply equal to the 
number of chiral edge modes. In the presence of inter- 
actions, for example in case of fractional Hall effect, the 
Green's function eigenvalues are no longer of this simple 
form. Instead, they generally have not only poles but 
also zeros as a function of iu, and No{k) counts the dif- 
ference of the signs of the poles and zeros (see Ref. [12] 
for details). 

Let us proceed to calculate the topological invariant 
Eq. (1) for a variety of fractional Hall states. As a 
warm up, let us first consider a simple Laughlin frac- 
tional Hall state corresponding to the Hall conductance 
CT^y = 1/ (2to -I- 1) and described by a ii'-matrix which re- 
duces to just one number, K = 2m + l. The edge Green's 
function in the position space for such a state is [13] 



G{x,t) 



1 



{x — vt) 



2m+l ' 



(8) 



where v is the velocity of its chiral excitations. Its Fourier 
transform is given by [13] 



G{n,k) = J 



^iQt—ikx 



dxdt 



{n + vkf 



(9) 



This can be obtained, for example, by introducing new 
variables = x ±vt in the integral. 
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We substitute the imaginary frequency w into Q, = iw 
to find 



G{w, k) 



{iuj + vk) 



(10) 



iLO — vk 

Then we plug this into Eq. (3). We can now take advan- 
tage of a simple relation, that for any function 

g = {iuj - Ay / {iw - Bf (11) 

with some real A and B, 



f 

J — ( 



^5 ^d^g = ^{13 s\g-n.B - asign A) (12) 



[which follows from Eq. (6)]. Then, putting Eq. (10) in 
Eq. (3) gives 



iVo(A) = -(2m+l)signA, (13) 



so that 



iV2 = 2m + 1 



(14) 



from Eq. (2). This is precisely the difference between the 
number of chiral modes and chiral zeros, represented by 
the denominator and the numerator of Eq. (10). 

First of all, N2 is indeed equal to the trace of iT- matrix, 
which in this case is just an odd integer. Also, we see 
that N2 is by no means equal to axy At the same time, 
this invariant is identical to the one computed for inte- 
ger quantum Hall effect with 2m -|- 1 filled Landau levels. 
Therefore, this invariant is not a unique identifier of a 
state, and does not necessarily change if a system under- 
goes a phase transition from an integer to a fractional 
Hall state. 

Let us now consider more general fractional Hall sys- 
tems described by a ii'-matrix. Their edge theory is given 
by the action [13] 

^ J [Kabdt<t>adxh - Vabdx(l)adxM ■ 

ah 

(15) 

Here Vab is the positive definite matrix of velocities of 
edge excitations, while K is the iiT-matrix, the matrix 
with integer entries which defines the topological order 
of the fractional Hall state. In this representation, we 
assume the operators creating fermions at the edge are 



(16) 



Let us compute the fermionic Green's function, substi- 
tute it into Eq. (2) and compute N2. 

To do that, it is advantageous first to note that since 
y is a positive definite symmetric matrix, we can pa- 
rameterize it in terms of some other symmetric matrix 



v = o\ 



(17) 



Changing the variables from (f> to (p = 0(f> (which now 
has units of \/ velocity) , we find the new action 

S=^ J dxdt [dtip Q-'^KQ-^dx^ - dx^dx^] ■ (18) 

Here the matrix notation for matrix products is used, for 
brevity. It is now convenient to diagonalize the symmet- 
ric matrix 



(19) 



where H is a diagonal matrix, and U is an orthogonal 
matrix, U^U = 1. Another change of variables 



brings the action to the simple form 
1 



(20) 



^^^^j dx<fiadt<fia/Va- {dx<Paf ■ (21) 



Here Va are the inverses of the diagonal entries of the 
diagonal matrix H, which obviously have units of velocity. 
It is straightforward to calculate the correlation functions 
of (fi now. They are given by 



{iPa{x,t) ^6(0,0)) 



^ , f diodk 6''=^-*"* 

-2Tndab / 77r-T2~n 

J (27r)^ cok/Va - 



k^' 



(22) 
(23) 



Doing the integral results in the expression 

X - Vat 

where L is the system size. 

Now we are interested in calculating 

Gab = (Mx,t) ^1(0,0)) = 

as follows from Eq. (16). To compute that, we work out 
(p in terms of ip to find 

J2KacMx,t)^ (^KmMO,0)^'^ = 
L 



OU' inlln 



x-t/n, 



uo 



(25) 



ab 



where the expression in the square brackets is understood 

as a product of matrices, and |n| In {L/ {x — t/W)) is un- 
derstood as a diagonal matrix with the diagonal entries 
\n{L/ [x — Vat)) / \va\- Finally, using the standard for- 
mula of Gaussian integration 



e 2 



i[<^<.(x,t)-</.6(0,0)]^ 



(26) 
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we find 



Sah 



(27) 



A word must be said about why tlie off-diagonal terms in 
Gab are zero. Evaluating these terms via Eq. (26) gives, 
for their L-dependence, 

Gab EcE<i('^<^<i°<io-^=<*^<*i')]^T^ . (28) 

If a ^ 6 the exponent is negative (columns of UO are 
linearly independent), so this goes to zero as L goes to 
infinity, justiiying the Kroncckcr delta in Eq. (27). If 
a = b, this is equal to 1 and Gaa does not depend on L 
as indicated in Eq. (27). Eq. (27) generalizes the Green's 
functions of fermions given by Eq. (8) to the case of a 
generic i^-matrix. 

We now need to perform a Fourier transform Gab{x, t). 
This is hard to do in general, and we expect that the an- 
swer is not as simple as it was in the previously consid- 
ered case Eq. (9). However, all we are interested in is a 
change in the argument of Goo(w, k) as uj is varied from 
—oo to 00, as explained in Eq. (6). This can be found in 
a relatively straightforward way. 

Let us calculate the Fourier transform of Eq. (27) by 
introducing the imaginary time t = it directly. 



Gaa{'^,k) = j drdx 



oioJT—ikx 



He {X + WXctY 



Here = {J^d UcdOdaf > 0, > is some positive 
constant with the units of velocity, and Xa = Va/v are 
some numbers of an arbitrary sign. Wc arc then inter- 
ested in calculating the change in the argument of G as 
k is kept fixed and uj is swept from — oo to oo. Let us 
show that this change is equal to 



N^{k) 



1 

2^ 



argGaa(a;,A;)||; 



-2-signfc, (30) 



where N^lk) is the contribution to No{k) from Gaa [it 
needs to be summed over a to find No{k)] and nia is 
defined as 

ma = Y,< signage = [OU^UUO] = Kaa- (31) 



To show that this is a diagonal entry of i^, we took ad- 
vantage of Eq. (19). nria is an integer since K is integer 
valued. 

Now we go back to evaluating Eq. (29). We change the 
variables from r, x to 



X + ivT = re ^. 



(32) 



This gives (for brevity, we suppress the index a in the 
Green's function and as an upper index of a) 



Jo vr^J_^ "^ric 



^iqr cos(0-0o) 



(cos(0) + iXcSm{(j))) 



— , (33) 



Here (pQ is the angle between the vectors {—k,Lu/v) and 
the vector (1, 0), r = Vx^ + vh^, q = ^/W+uPjv^, and 
n = J2a ^a- rewrite this integral as 



G 



rdr 



vr" 



Further analysis depends on whether Xc are positive or 
negative. If Xc is positive, we take e**^'^+'°"'(l+.-rc)/2 out- 
side the appropriate bracket, otherwise we do the same 
with e-*('^+'^»)(l - Xc)/2. We find 



G' 



-im(po 



r rdr r 
J vr" J_^ 



d(f>- 



rjiqr cos{(f))—im4 



+ g-2i((^+(?io)signxcp^^ ' 

(35) 

where m was defined above in Eq. (31) and = (1 — 
|a:,|)/(l + |a:,|), \p,\<l. 

As uj is swept from — oo to oo, i/io grows from — 7r/2 
to 7r/2 if A; < and decreases from 37r/2 to tt/2 ii k > 
0. We would like to know the change in the argument 
of G during this process. The outer explicit factor of 
g-im0o obviously contributes Trmsignfc to that change. 
Let us argue that the rest of the expression does not 
contribute at all to the change of the argument. To do 
that, we rewrite everything which multiplies e' 
this expression as 



-imcpQ 



f rdr r 



piqr cos{4>)—im4' 



_|_ g-2i((^+<^o) sign Xc ^ 



(36) 



It is clear that if all = (which corresponds to all 

velocities being equal to each other in magnitude and 
equal to v), then this expression is independent of (po 
and its argument does not wind at all as a function of 
(^0- As pc are increased, in order for the argument to start 
winding, one should be able to find such (/iq, at some pc, 
that this expression is equal to zero (or infinity, but this 
expression is finite as long as all |pc| < 1- 

However, if this expression is zero, that means that the 
Green's function has a zero at those values of pc- Zeros 
of the Green's function can occur only at cj = (see page 
168 of Ref. [14] or Ref. [12] for a slightly stronger version 
of this statement with a more detailed discussion), a; = 
corresponds to 0o = or 00 = it. However, if such a zero 
were to occur, then for a related system with velocities 
v'a = v/va where p'^ = —pc, the zeros would occur at 
= 7r/2 which contradicts the theorem that zeros occur 
only at w = or w = TT. Therefore, the end result is that 
this is impossible. We conclude from this that Eq. (36) 
does not wind around zero of the complex plane as (po 
changes from —n/2 to tt/2, or from 37r/2 to tt/2. 

Taken together, this shows that indeed the contribu- 
tion to No{k) from Gaa is given by Eqs. (30) and (31) 
to be Kaa- Summing over all the entries of the Green's 
function we find that 



Noik) 



1 



; sign k ■ tiK. 



(37) 
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That, in turn, leads to the anticipated answer, 
N2 = No{A)-No{-A)=tvK. 



(38) 



Note that this expression is not invariant under the 
change of basis (going from the basic fields (j)a in Eq. (15) 
to some linear combination thereof). Indeed, the topo- 
logical invariant is defined in the preferred basis where 
the fermion operators are simply defined, as in Eq. (16). 
It is possible to rewrite the expression for the topological 
invariant in the basis independent form with the help of 
the "electron lattice" matrix C introduced in Ref. [15]; 
since, in the present context, the purpose of this matrix 
would just be to take the matrix K back to the preferred 
"electron" basis, we omit this. 

Now, let us now briefly look at the non-Abelian Read- 
Rezayi states [10] (which include the Pfaffian, or Moore- 
Read, state as a particular case [11]). The edge theory of 
these states include one charged and one neutral particle 
(a Luttinger liquid-type boson and a parafermion, which 
reduces to a Majorana fermion in case of the Moore- Read 
state). The Green's functions at the boundary can be 
easily derived from conformal fleld theory and are equal 
to 



1 



{x - Vnt) ™ (x - Vet) 



M-l-i 



(39) 



Here N is the level of the Read-Rezayi state {N = 2 
corresponds to the Pfaffian state) and M is an odd posi- 
tive integer (even M would correspond to a bosonic state 
which we do not discuss here). The filling fraction of 
these states is known to be [10] 



N 



MN + 2 



(40) 



(this number represents the filling fraction of the partially 
filled Landau level; in the presence of lower completely 
filled Landau levels the filling fraction can be larger than 
this number by an integer). 

It is now straightforward to use Eqs. (29) and (30) 
to find that the topological invariant is the sum of two 
exponents in Eq. (39), or 



N2 = M + 2. 



(41) 



Interestingly, it is independent of N. For the 5/2 plateau 
in fractional Hall effect, understood as M = 1, A'' = 2 
Pfaffian state, this gives A^2 = 3, the same as for the 
Laughlin = 1/3 state. 

A final non-Abelian state of interest to us is the anti- 
Paffian state ai u = 5/2 [16-18]. Its edge Green's func- 
tion is given by 



G 



1 



{X - Vet) {X + Vnt) 



(42) 



Importantly, here the charged and neutral modes are 
counterpropagating. As a result, employing Eqs. (29) 



and (30) we find that the topological invariant is the dif- 
ference of the two exponents in Eq. (42) or 



No = 1, 



(43) 



the same as in the simple axy = e'^/h integer quantum 
Hall state. Note that in all these cases, the invariant 
N2 appears simply to reflect the scaling dimension of the 
electron operator [13]. 

While the value of the invariant would not be helpful 
in trying to distinguish the anti-Pfaffian from a simple 
integer Hall state, it can scirvc to distinguish the Pfaffian 
(Ar2 = 3) from the anti-Pfaflaan (ATa = 1) 1/ = 5/2 state. 

In conclusion, we have defined an invariant N2 for frac- 
tional Hall states via Eq. (1) and calculated it for a va- 
riety of fractional and non-Abelian quantum Hall states. 
We find that the invariant is not a unique identifier of the 
state; a phase transition can occur while the invariant 
may stay the same. However, within the standard low- 
energy descriptions employed here, it appears that the 
invariant cannot change without a phase transition since 
its value is defined by the state; this should be contrasted 
with the situation in one dimensional systems where the 
invariant can change without a phase transition [19]. 

We would like to add one final remark: it is possible to 
try to define the invariant not via the Green's functions 
but rather via the phases (boundary conditions) across 
the system [4] which can lead to an expression equal to 
the fractional Hall conductance, unlike the invariant con- 
structed here, which is always an integer and is not equal 
to the Hall conductance. This may be an interesting av- 
enue to pursue in the future. 

VG would like to acknowledge NSF grant no. DMR- 
1205303 for support, and is grateful to A. Altland for 
useful discussions, in particular for pointing out the di- 
agonal structure of the Green's functions in Eq. (27). 
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